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INTRODUCTION 


A  flexible-fin  decelerator  to  orient  and  stabilize  "smart"  submunitions  has  been  developed 
and  patented  (ref  1).  These  submunitions  are  ejected  from  a  rapidly  spinning  projectile  and 
descend  vertically  over  a  battlefield,  searching  for  armored  targets  (fig.  1).  The  flexible  fin 
drives  the  submunition  in  a  lunar  motion  at  a  predetermined  steady  state  spin  rate  and  fall 
velocity.  A  lunar  motion  is  required  for  target  scan  and  for  the  sensor  to  lead  the  warhead  axis. 

Since  the  submunitions  are  ejected  from  a  rapidly  spinning  projectile  (up  to  260  Hz),  reliable 
deployment  of  the  flexible  fin  could  be  a  problem. 

The  objective  of  the  present  work  is  to  develop  a  numerical  simulation  to  model  the  deploy¬ 
ment  of  a  flexible-fin  decelerator  from  a  rapidly  spinning  body.  Such  a  simulation  allows  for  the 
evaluation  of  various  fin  packing  techniques  and  fin  strength  without  resorting  to  the  usual 
cut-and-try  test  procedures. 

The  method  uses  an  approach  similar  to  that  of  Purvis  (refs  2  and  3)  in  which  the 
decelerator  is  modeled  as  a  series  of  masses  connected  by  flexible  elements.  The  motion  of 
the  masses  is  determined  solely  by  tension  forces;  aerodynamic  forces  will  be  included  in  future 
efforts.  Further  limitations  to  this  study  are  constant  submunition  spin  rate  and  two-dimensional 
motion  confined  to  the  plane  of  the  top  of  the  submunition. 


DISCUSSION 

The  system  to  be  simulated  consists  of  the  submunition  body,  the  flexible  fin,  and  the  tip 
weight  (fig.  2).  The  fin  is  modeled  as  a  series  of  masses  connected  to  each  other  (and  to  the 
body  and  tip  weight)  by  springs  and  dampers.  A  drawing  of  the  system  model  with  the  fin 
stretched  out  is  shown  in  figure  3.  The  undeflected  length  of  the  fin  elements  is  s  . 

Of 


The  equations  of  motion  of 

a  given  mass  are: 

m'x.  =  C 

T  +  D.Jcose.  -  c 

fj  ,  +  Di.1)cos0.  , 

(1) 

my.  =  C 

F.  +  D.)sin9.  -  (7 

'i-i  +  Di-i)sin0i-i) 

(2) 

The  coordinate  system  and  angular  definitions  are  shown  in  figure  4.  The  only  forces  assumed 
to  act  on  the  masses  are  tension  and  damping  forces.  Tension  and  damping  between  masses  i 
and  i-1  are  T.  and  Djt  respectively.  Tension  is  computed  by  calculating  the  distance  between 

the  masses  (s^  and  comparing  it  to  the  undeflected  distance  (so).  If  Sj  is  less  than  so,  tension  is 

zero.  The  spring  constant  is  calculated  using  the  breaking  strength  and  percent  elongation  at 
breaking  of  the  fin  material: 


Tension  is  given  by: 


T  -  (s,  -  s0)k  (4) 

Damping  is  assumed  linearly  proportional  to  the  velocity  of  one  mass  relative  to  the  other. 
The  distance  between  masses  is  given  by: 

Si  =  V(Xi  -  Xj.,)2  +  (yj  -yj-i)2  (5) 


The  derivative  of  this  quantity  is  the  relative  velocity: 

=  l(xi  -Xj.QtXi  -  Xi - 1 )  +  (yi  -yi-i)(yi  -yi-i)l 

V(Xi  -Xi.,)2  +  (yi  -  yj  - 1  )2 


(6) 


Therefore,  damping  is  given  by: 


Di  =  cdsi 


(7) 


where  Cd  is  the  damping  coefficient.  As  with  the  tension,  if  Sj  is  less  than  so,  damping  is  zero. 


A  flow  chart  of  the  deployment  simulation  is  shown  in  figure  5.  A  fourth  order  Runge-Kutta 
routine  is  used  to  solve  the  equations  of  motion  for  each  mass. 


RESULTS 

Several  test  cases  were  run  to  demonstrate  the  simulation  and  evaluate  three  different  fin 
packing  techniques.  The  submunition  weight  was  10  lb  and  had  a  diameter  of  5  in.  The  spin 
rate  chosen  was  450  rad/s,  representative  of  the  spin  rate  at  cargo  ejection  of  a  low  zone 
155-mm  firing.  The  fin  had  a  10-in.  span,  3-in.  chord,  weighed  0.02  lb,  and  was  divided  into  4 
masses  (0.005  lb  each).  The  tip  mass  weighed  0.2  lb.  For  these  preliminary  cases,  the  damp¬ 
ing  coefficient  (Cd)  was  set  to  zero.  The  fin  breaking  strength  was  assumed  to  be  800  lb  and 

elongation  at  breaking  (eBr)  25%. 

When  dealing  with  relatively  small  masses  and  large  spring  constants,  system  frequencies 
can  be  high.  Choosing  a  small  enough  timestep  is  crucial  if  large  errors  are  to  be  avoided.  The 
smallest  mass  and  largest  spring  constant  can  be  used  to  calculate  the  highest  frequency 
possible  in  the  system.  This  frequency  is  given  by: 
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The  maximum  spring  constant  is  obtained  from  equation  3  using  the  undeflected  length  be¬ 
tween  fin  masses  (so),  in  this  case  2.0  inches.  From  equation  3,  kmax  =  19,200  Ib/ft.  The 

smallest  mass  in  the  system  is  a  fin  mass,  0.005  lb.  From  equation  8,  w  =  1 1 ,120  rad/s.  The 
period  is  0.00056  sec. 

To  test  the  effect  of  the  size  of  the  timestep,  the  simple  case  of  two  fin  masses  connected 
with  a  2.0-in.  long  spring  was  simulated.  At  time  zero,  one  of  the  masses  is  given  an  initial 
velocity  of  94  ft/s.  The  position  and  velocity  versus  time  of  this  mass  can  then  be  compared  to 
the  analytic  solution.  The  excat  solution  is  compared  to  the  simulation  for  the  indicated  time 
step  in  table  1 .  Since  the  velocity  varies  sinusoidally  about  47.0  ft/s,  the  difference  in  velocity  is 
normalized  with  the  average  velocity,  47.0  ft/s.  The  difference  in  position  is  normalized  with  the 
exact  position.  For  a  timestep  of  0.000025  sec,  the  simulated  position  and  velocity  deviates 
less  than  2%  from  the  exact  solution.  For  the  test  cases  described  below,  a  timestep  of 
0.000025  sec  was  used,  roughly  20  integrations  per  cycle. 

The  first  simulated  packing  technique  was  lumping  an  of  the  fin  masses  and  tip  weight  at 
the  fin  root  (fig.  6a).  The  tip  weight  was  given  the  tangential  velocity  at  the  outside  edge  of  the 
submunition  (S/M)  and  was  free  to  deploy.  When  the  tip  mass  moved  far  enough  away  from  the 
S/M  to  put  tension  between  it  and  the  first  fin  mass,  the  first  fin  mass  was  allowed  to  deploy. 
The  remaining  fin  masses  were  allowed  to  deploy  sequentially  in  the  same  manner.  The 
position  of  each  mass  and  the  submunition  was  plotted  with  a  graphics  program.  A  sequence  of 
these  plots,  each  at  a  later  time,  is  shown  in  figure  7.  The  deployment  sequence  shows  the  fin 
falling  behind  the  submunition  and  starting  to  cross  over  the  top.  In  a  tactical  environment, 
aerodynamic  loads  on  the  fin  (not  included  in  the  present  simulation)  as  it  is  deploying  may  tend 
to  allow  it  to  "catch  up"  to  the  submunition.  Tension  in  the  flexible  members  between  the 
masses  is  shown  in  figure  8.  The  maximum  tension  calculated  was  about  1250  lb  and  occured 
between  the  second  and  third  fin  masses.  For  this  packing  technique,  it  is  likely  that  the  fin 
would  fail  at  deployment. 

The  second  simulated  packing  technique  was  to  position  the  fin,  stretched  out,  around  the 
periphery  of  the  submunition  wrapped  in  the  spin  direction  (fig.  6b).  All  of  the  masses  were 
given  the  tangential  velocity  at  the  outside  edge  of  the  submunition  and  were  allowed  to  deploy 
at  time  zero.  The  resulting  deployment  sequence  (fig.  9)  shows  the  fin  to  deploy  "ahead"  of  the 
submunition  initially  and  become  fully  erect  at  about  0.012  seconds.  The  fin  then  falls  behind 
the  submunition,  but  in  a  tactical  environment  and  once  its  erect,  aerodynamic  forces  will  keep  it 
erect.  Tensions  are  shown  in  figure  10.  Maximum  tension  is  about  750  lb  and  occurs  between 
fin  mass  four  and  five.  The  fin  strength  is  adequate  for  this  packing  technique. 

The  third  simulated  packing  technique  was  to  position  the  fin,  stretched  out,  around  the 
periphery  of  the  submunition,  wrapped  opposite  to  the  spin  direction  (fig.  6c).  All  of  the  masses 
were  given  the  tangential  velocity  at  the  outside  edge  of  the  submunition  and  were  allowed  to 
deploy  at  time  zero.  The  resulting  deployment  sequence  (fig.  1 1 )  shows  that  the  fin  immediately 
crosses  the  top  of  the  submunition,  becomes  partially  erect,  and  then  falls  behind.  It  is  likely 
that  the  deployment  in  this  manner  will  result  in  a  twisted  fin.  Tensions  are  shown  in  figure  12. 
Maximum  tension  is  about  1400  lb  and  occurs  between  mass  five  and  the  submunition.  Be¬ 
sides  twisting,  it  is  likely  that  the  fin  would  fail  at  deployment. 
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The  packing  technique  resulting  in  the  cleanest  and  lowest  deployment  loads  is  wrapping 
the  fin  around  the  submunition  periphery  in  the  spin  direction. 


CONCLUSIONS 

A  simulation  to  model  the  deployment  of  a  flexible  fin  decelerator  has  been  developed.  At 
present,  motion  is  confined  to  two  dimensions  and  determined  solely  by  tension  forces.  An 
obvious  extension  to  this  work  is  to  add  aerodynamic  forces  and  to  allow  motion  in  three  dimen¬ 
sions.  Also,  the  simulation  must  be  validated  by  test  data. 

Of  the  test  cases  simulated,  the  most  promising  deployment  resulted  when  the  fin  was 
wrapped  around  the  periphery  of  the  submunition  in  the  same  direction  as  the  submunition  spin. 
This  resulted  in  a  clean  deployment  with  the  lowest  fin  tension  forces. 


Table  1.  Comparison  of  simulation  to  exact  solution 


At  =  0.0001  At  =  0.00005 _  At  =  0.000025 


t  (sec) 

X 

exact 

(ft) 

• 

X 

exact 

(ft/s) 

x  . 

sim 

error 

(%) 

*sim 

error 

(%) 

x  . 

sim 

error 

(%) 

- f - 

X  . 
sim 

error 

(%) 

X 

sim 

error 

.  (%)._ 

x  . 

xim 

error 

(%) 

0.0001 

0.007689 

46.88 

0.007463 

2.94 

47.86 

-2.1 

0.007680 

0.12 

47.10 

-0.47 

0.007688 

0.013 

46.93 

0.11 

0.0002 

0.009385 

0.0006 

0.009501 

-1.24 

6.85 

-14.6 

0.009413 

-0.30 

0.28 

-0.59 

0.009391 

-0.06 

0.01 

-0.02 

0.0003 

0.11111 

47.36 

0.011742 

-5.68 

44.79 

5.5 

0.011138 

-0.24 

47.31 

0.11 

0.011112 

-0.009 

47.21 

0.32 

0.0004 

0.018831 

94.00 

0.018627 

1.08 

81.24 

27.1 

0.018774 

0.30 

93.43 

1.2 

0.018818 

0.07 

93.98 

0.04 

0.0005 

0.026489 

46.40 

0.025510 

3.70 

50.14 

-8.0 

0.026444 

0.17 

47.51 

-2.4 

0  026487 
0.008 

46.65 

-0.53 

0.0006 

0.028154 

0.0056 

0.028422 

-0.95 

17.84 

-37.9 

0.028238 

-0.30 

0.85 

-1.8 

0.028174 

-0.07 

-.03 

-0.05 

0.0007 

0.029912 

47.84 

0.031190 

-4.27 

43.24 

9.8 

0.029974 

-0.21 

46.30 

3.2 

0.029914 

-0.007 

47.48 

0.77 

0.0008 

0.037661 

93.99 

0.037348 

0.83 

71.79 

47.2 

0.037549 

0.29 

92.87 

2.3 

0.037635 

0.09 

93.96 

0.06 

0.0009 

0.045288 

45.91 

0.043752 

3.40 

51.12 

-11.1 

0.045208 

0.18 

47.90 

-4.2 

0.045286 

0.004 

46.38 

-1.0 

0.0010 

0.046923 

0.015 

0.047269 

-0.73 

25.96 

-55.2 

0.047063 

-0.30 

1.41 

-3.0 

0.046956 

-0.07 

0.05 

-0.07 

0.0011 

0.048712 

48.33 

0.050468 

-3.60 

42.71 

12.0 

0.048810 

-0.20 

45.91 

5.1 

0.048715 

-0.006 

47.76 

1.2 

0.0012 

0.056492 

93.98 

0.056125 

0.65 

64.83 

62.0 

0.056324 

0.30 

92.30 

3.6 

0.056453 

0.07 

93.94 

0.09 

.'«V< 


vm 

*.*V 


m 


iS 


-CALCULATE  FORCES 
-INTEGRATE  FOR  NEW 
POSITIONS  &  VELOCITIES 


-COMPUTE  TENSION  AND 
CHECK  IF  TIME  STEP 
NEEDS  ADJUSTMENT 


-TEST  FOR  NEXT  MASS 
DEPLOYMENT 
(IF  NECESSARY) 


-WRITE  OUTPUT 


Figure  5.  Simulation  flowchart 


Tensions  versus  time  for  sequential  deployment  from  root 


Hgure 


t  =  0.008S 


=  0.010s 


Figure  1 1 .  Deployment  sequence,  fin  wrapped  opposite  to  spin  direction 
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Figure  11.  (cont) 


Tensions  versus  time  for  fin  wrapped  opposite  to  spin  direction 


Figure  12.  (cont) 
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Damping  force  between  masses  i  and  i.  lb 
Breaking  strength  of  fin  material,  lb 
Maximum  spring  constant  in  system,  Ib/ft 

Spring  constant,  Ib/ft 
Mass  of  element  i,  lb 

Minimum  mass  in  the  system,  lb 
Distance  between  mass  i  and  i  +  1 ,  ft 
Unstretched  distance  between  mass  i  and  i  +  1 ,  ft 
Tension  between  mass  i  and  i  +  1 ,  lb 
Coordinates  of  mass  i,  ft 

Angle  (with  respect  to  horizontal)  of  line  segment  connection  mass  i  and  i  +  1 , 
deg 

Maximum  frequency  in  system,  rad/s 
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CLOSED  FORM  SOLUTION  OF  SIMPLE  TWO  BODY  PROBLEM 
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Two  bodies  of  mass  mt  and  m2  are  connected  by  a  massless  spring  of  free  length  X  and 

spring  constant  K  (fig  A1).  Free  body  diagrams  of  the  masses  are  also  shown.  The  equations 
of  motion  of  the  two  bodies  (assuming  no  damping)  are: 

"V,  =  k(x2  -  x,  -£)  (A1 ) 

m2x2  = -k(x2  -  x,  -%)  (A2) 

These  equations  can  be  rewritten  as: 

x,  +  w/x,  =  wt2(x2  -l)  (A3) 

”  .  ...  2 . 2,  .  n  .  .... 


x2  +  w22x2  =  w22(x2  +  l ) 


where 


wi,2  V  m,.2 


Introduce  a  differential  operator  P  into  equation  (A1): 


P2x,  +  w2x,  =w12(x2-/t) 


solve  for  x  ■ 


W?<X2  -JO 


X'  ”  (P2  +  wf  ) 


Substitute  equation  (A6)  into  equation  (A4): 


X2  +  w2  x2  = 


2rw?(x2-ju  0i 


Equation  (A7j  can  be  simplified  to: 
'x2  +  (wt2  +  w22)x2  =  0 


Define  y  =  x2  and  w2  =  w1 2  +  w22  so 


y  +  w  y  =  0 


The  solution  of  equation  (A9)  is  of  the  form 
y  =  A  sin  wt  +  B  cos  wt 
Since  y  =  x 


x2  =  A  sin  wt  +  B  cos  wt 


Integating  once 


x2  =  -A/w  cos  wt  +  B/w  sin  wt  +  C 


Integrating  again 


x2  =  -A/w2  sin  wt  -  B/w2  cos  wt  +  ct  +  D 


(A10) 


(AH) 


(A12) 


(A13) 


Initial  conditions  may  now  be  used  to  determine  the  constants  A,  B,  C,  and  D.  For  the  test 
case,  assume  x^O)  =  0,  x2(0)  =  x^O)  =  0,  and  x2(0)  =  Vq.  Substituting  for  x^O)  and  x2(0)  in 

equation  (A2),  x2(0)  =  0.  If  x‘2(0)  =  0,  from  equation  (All),  B  =  0.  So  far  the  solution  is: 


x2  =  -A/w2  sin  wt  +  ct  +  D 


Since  x2(0)  =  from  (A  14),  D  =  J^.  If 


x2  =  A  sin  wt 


(A14) 


(A15) 


then 


x2  =  Aw  cos  wt 


Form  equation  (A2), 

x'2  = -w22(x2  -  Xi  -JO 


(A16) 


(A17) 


so 


X2  =  W2  (X2  ‘  Xl) 


(At  8) 


Substituting  for  xt(0)  and  x2(0) 


Finally,  substituting  x2(0)  into  equation  (A12) 


V 

0 


V0  +  c 


so 


w; 

c  =  V  (1  -  -§) 

O'  7 


The  solutuion  of  equation  (A8)  is: 

o  2 

Wo  Wo  v  A 

x,=  4v  sin  wt  +  V  (1  -  t  +  X 

2  yy3  o  o'  yy«:  '  ^ 


and 


.  w; 

2  w2 


w; 


vocoswt  + V0(1  -  ^5) 


(A20) 


(A21) 


(A22) 


(A23) 
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